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ABSTRACT 

Low-frequency Alfven-wave turbulence causes ion trajectories to become chaotic, or "stochastic," when the 
turbulence amplitude is sufficiently large. Stochastic orbits enable ions to absorb energy from the turbulence, 
increasing the perpendicular ion temperature T±i even when the fluctuation frequencies are too small for a 
cyclotron resonance to occur In this paper, an analytic expression for the stochastic heating rate is used in 
conjunction with an observationally constrained turbulence model to obtain an analytic formula for T±i as a 
function of heliocentric distance r, ion mass, and ion charge in coronal holes at 2Rq < r < 157^0. The resulting 
temperature profiles provide a good fit to observations of protons and 0+^ ions at 2Rq < r < 3Ra, from the 
Ultraviolet Coronagraph Spectrometer (UVCS). Stochastic heating also offers a natural explanation for several 
detailed features of the UVCS observations, including the preferential and anisotropic heating of minor ions, 
the rapid radial increase in the 0+^ temperature between 1 .6/?0 and 1 SRq,, and the abrupt flattening of the 0+^ 
temperature profile as r increases above \ .9Rq. 

Subject headings: solar wind — Sun: corona — turbulence — waves — MHD 



1. INTRODUCTION 

Measurements of ion and electron temperatures show 
that the solar wind undergoes substantial heating through- 
out its transit from the solar surface to the heliospheric 
termination shock. A number of considerations suggest 
that much of this heating results from Alfven- wave tur- 
bulence, as suggested over forty years ago by 'C olemanI 
(|T968). For example, Alfven-wave-like fluctuations in 
the velocity, magnetic field, and electric field are often 
present in the interplanetary medium, with power spec- 
tra that vary approx imately as power laws over a broad 
range of len gth scales ("Belcher & Davisll971J: lGoldste in etaT 



1995 



lTu & Marsch 1995; Bruno & Carbonel2 005: Ba le et al 
2005|). Alfven waves hav e also been measured re- 
motely in the chrom osphere jDe Pontieu et al.l l2007h and 
corona (Tomcz vk et alj|2007l) . and radio observations at he- 
liocentric distances r as small as ^ 5Rq, reveal a broad 
spectrum of density fluctuations that are consistent with 
passive-scal ar mixing of entropy modes by Alfven-wav e 
turbulence jHarmon & ColesI 120051: iChandran et"al] 120091) . 
At r > 6O/?0, the amplitudes of velocity and magnetic- 
field fluctuation s are positively correlated with tempera- 
ture dGrappin et al. 1990), and the turbulent dissipation rate 
implied by the measured fluctuation amplitudes is compa- 
rable to the heating rate inferred from m easurements of the 
proton and electron temperature profiles jSmith et alj|2001l : 
ICranmer et al. 2009; Stawarz et al. 2009). 

On the other hand, it has for some time been unclear 
whether Alfven-wave turbulence can explain the ion temper- 
atures and temperature anisotropics observed in the fast so- 
lar wind and coronal holes (the open-magnetic-field-line re- 
gions from which the fast wind emanates). For example, 
observations from the Ultraviolet Coronagraph Spectrome- 
ter (UVCS) show that minor ions such as 0+^ are heated to 
temperatures greatly exceeding the proton and electron tem- 
peratures at 2Rq < r < 3/?0, with thermal motions that are 
much more rapid in directions perpendicular to the magnetic 
field B than along the magnetic field — i.e., 7x ^ Tii (Kohl 



et al. 1998; Li et al. 1998; Antonucci et al. 2000). Simi- 
larly, in-situ measurements from the Helios and Wind satel- 
lites show that 7xp > 7]|p > T^. in low-p fast-wind streams at 
6Q)Rq, < r < 1 AU, where Txp and 7||p are the perpendicular 
and parallel proton temperatures, Tg is the electron tempera- 
ture, and P = Snp/B^ is the ratio of the plasma pressu re p to 
the m agnetic pressure (Phillips & Gosling 1990; Marsc het alj 
I2OOI . 

These observations lead to an apparent difficulty for solar- 
wind-heating models based on Alfven-wave turbulence for 
the following reasons. In Alfven-wave turbulence, the energy 
cascade is anisotropic, transporting energy primarily to small 
scales measured perpendicular to B rather than small scales 
along B (Shebalin etal. 1983; Goldreic h & Sridhar 19951; 
Ng & Bhattacharied 1 1 9961; iGaltier et airifoOoTlCho et al J 



2002ir In wavenumber space, energy cascades primarily to 
larger k±, and only weakly to larger ^||, where k± and k^^ are 
the components of the wave vector k perpendicular and par- 
allel to B. At perpendicular scales X± of order the proton gy- 
roradius pp, the Alfven-wave c ascade transitions to a kinetic- 
Alfven -wave (KAW) casc ade (iBale et al1l2005l: iHowes et alj 
2008b; iSchekochihineFal . 2009; Sahra ouiet all 120091) . and 
at X± < pp the fluctuations dissipate. The rms amplitude of 
the magnetic-field fluctuations at X± ^ Pp is <C B, and thus it 
seems reasonable to assume that the KAW fluctuations damp 
at the same rate as linear KAWs at the same k. For the p val- 
ues found in coronal holes and the solar win d, linear KAWs 
undergo significant e lectron Landau damping ('Ou ataerll 19981 ; 
iLeamon et al.lll99"9l) . However, the KAWs produced by the 
anisotropic cascade do not undergo ion cyclotron damping, 
because their frequenci es (0 are much less than the pro- 
ton cyclotron frequ ency (ICranmer & van Ballegooiienll2003l ; 
IHowes et alll2008ij) . In addition, when P < 1, flie ion ther- 
mal speeds are <C va, where va is the Alfven speed. Ions are 
thus unable to satisfy the resonance condition co — fe||V|| = 
for Landau damping or transit-time damping, where V|| is 
the velocity component parallel to B, because ft)/A;|| > va for 
KAWs ('O uataertlfTggI iHoUwegl [19991) . As a consequence, 
linear damping of KAWs by ions is neghgible when p <C 1, 
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suggesting that low-frequency Alfven-wave/KAW turbulence 
results in negligible ion heating in coronal holes and low-p 
fast- wind streams. 

A number of studies have gone beyond the framework of 
linear wave theory to explore the possibility of perpendicu- 
lar ion heating by low-frequency turbulence in the solar wind 
(e.g., Dmitruk et al. 2004; Markovskii et al. 2006; Bourouaine 
et al. 2008; Parashar et al. 2009; Lehe et al. 2009; - see 
also Drake et al. 2009). Several investigations have shown 
that low-frequency Alfven waves and KAWs can indeed 
cause perpendicular ion heating if the wave amplitudes are 
sufficiently large (Johnson & Cheng 2001; Chen et al.H200Tl: 
IWhite et al.l 12002; Voitenko & Goossens 1004). This type 
of heating is some times referred to as "stochastic heating," 
since it is predicated upon the stochastic ion orbits that arise 
in the presence of large-amplitude fluctuations in the elec- 
tric or magnetic field. This paper addresses the question of 
whether stochastic ion heating can explain the anisotropic ion 
temperatures that are observed in cor onal holes. The analy- 
sis is based on the recent studies by IChandran et al.l (120101) 
and|Chandran & Hollweg (120091) . which are reviewed in sec- 
tions|2]and[3] Stochastic heating in coronal holes and the inner 
solar wind is then discussed sectionS) 

2. STOCHASTIC ION HEATING 



kinetic energy in the absence of cooling: 

fh^er'nriexp(^^j, (2) 

where ci is a dimensionless constant. An important conse- 
quence of equation (|2]i is that stochastic heating is inherently 
self-limiting — as T\x grows, Z\ decreases, and the stoc hastic 
heating rate decreases as a result. IChandran et al.l (120101) eval- 
uated C2 by simulating test-particle protons interacting with a 
spectrum of randomly phased Alfven waves and KAWs in a 
low-p plasma and found that C2 = 0.34. However, they argued 
that C2 is smaller for Alfven-wave/KAW turbulence than for 
randomly phased waves, because turbulence p roduces coher- 
ent st ructures that increase orbit stochasticity (iDmitruk et alj 
120041) . In the present paper, C2 is taken to be the same for all 
ion species, and the value 

C2 = 0.15 (3) 

is chosen in order to match UVCS observations of 
temperatures, as discussed further in section |4] below. 
IChandran et al.l (12010 ) also showed that in low-p plasmas 
stochastic heating by low-frequency Alfven-wave/KAW tur- 
bulence primarily increases an ion species' perpendicular 
temperature Txi, rather than its parallel temperature Tli;. 



If the amplitudes of Alfven waves or KAWs with A-x ^ Pi 
and CO < i2i are sufficiently large, then ions undergo "stochas- 
tic heating," where A-x is the wavelength measured perpen- 
dicular to B, CO is the wave frequency, Pi ~ vxi/iii is the 

ion gyroradius, vxi ~ \J Ik-^ 7xi / m\ is the rms value of vx 
(the ion velocity component perpendicular to B in the local 
plasma frame), = q\B lm\c is the ion cyclotron frequency, 
and mi, q\ and T^ m are the ion mass, charge, and perpendic - 
ular temperature (Johnson & Cheng 20011; iChen et al.ll2001[) . 
IChandran et al. (2010) showed that the amplitude threshold 
for strong stochastic heating when p < 1 can be expressed in 
terms of the quantity 



8vi 

Vxi' 



(1) 



where 5vi is the rms amplitude of the velocity fluctua- 
tion at A-x ^ Pi. For Alfven-wave and KAW fluctuations. 



Ei - ^i5<I>i/;i 



"Xi' 



where 5<I>i is the rms amplitude of the 



electrostatic -potential fluctuation at Ax ^ Pi, and thus 8i is 
roughly the fractional change in an ion's kinetic energy in- 
duced by gyroscale fluctuations during the course of a single 
cyclotron period. If Ei is sufficiently small, then a thermal 
ion's orbit in the plane perpendicular to B closely approxi- 
mates a closed circle in some suitably chosen reference frame. 
In this case, the ion's magnetic moment /j = miv^ /2B is nearly 
conserved, and perpendicular ion heating is weak dKruskall 
119621) . On the other hand, as Ei increases towards unity, a ther- 
mal ion's orbit becomes increasingly chaotic, /j conservation 
breaks down, and perpendicular ion heating becomes increas- 
ingly strong. For protons, vxi ^ P'/^va, 5vi ^ vxhB\/B, and 
Ei - p-'/^SBi/B, where 5Bi is the rms amplitude of the mag- 
netic fluctuation at Ax ~ Pi. Thus, when p ^ 1, Ei is much 
greater than bB[ j B. 

Usin g phenomenological arguments, IChandran et al.l 
( 120101) derived an expression (their equation (37)) for the 
time scale ?h on which stochastic heating doubles an ion's 



3. ALFVEN-WAVE TURBULENCE IN CORONAL HOLES 

Convective motions at the surface of the Sun randomly stir 
the footpoints of open magnetic field lines, launching low- 
frequency Alfven waves that propagate into coronal holes and 
then on into the solar wind. Radial variations in the Alfven 
speed VA{r) couple the outward-propagating Alfven waves 
to inward-propa gating Alfven waves (Heinemann & Olber^ 
[1980; Vein 1993riHollweg & Isenberg|l2007HCran mer 2010), 
and nonlinear interactions among counter-propagating waves 
then cause the wave energy to cascade to small scales 
and dissipate, heating the ambient plasma (IVelli et al."1989t 
patthaeus et al. 1999; Cranmeretal. 2007; Verdini et aQ 
vOlO). The value of the correlation length or "outer scale" 
at the coronal base, denoted L±b, is in the present study taken 
to be 5 X 1 0^ km, com parable to the spacing of photospheric 
flux tubes ( Spruitl[T98 1). At larger r, the value of Lj_ is taken 
to increase in proportion to the radius of a flux tube, so that 



Boin,) 



Bo{r 



1/2 



(4) 



where Bq is the background magnetic field strength and r\, ~ 
IRq corresponds to the coronal base (i.e., x = 1 in equation (|9]l 

b elow). 

IChandran & H ollweg (2009) developed an analytical model 
of Alfven-wave reflection and turbulent heating based upon 
a phenomenological treatment of the energy cascade. Using 
this model, they calculated the radial profiles of the turbulent 
dissipation rate per unit mass F and rms wave amplitude at 
the outer scale , denoted 8vo, extending the previous model of 
IDmitruk et al.l ( 12002.) to account for the solar-wind velocity. 
For the case that mos t of the Alfven-wave en ergy is at periods 
of ~ 1 hour or longer. lChandran & Hollweg (12009) found that 
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Fig. 1. — The short-dashed line is the rms amplitude of the fluctuating 
velocity from equation (5j with Si'a = 155 km/s, where n, Bq, and U are taken 
from equations (sj through jlOt . The solid-line box represents the range of 
non-thermal line widths obtained by Esseretal. 1 1999). The dotted line is 
the solar-wind velocity from equatio njlOt . and the long-dashed hne is the 
Alfven speed from equations (8) and (9f 



and 



r = (i+ii-'/2)5v2 



dvA 



dr 



(6) 



for < r < 2QRq, where rm is the radius at which the Alfven 
speed v'A = Bq/ ^JAnnnip obtains its maximum value, n is the 
proton number density, nip is the proton mass, 5va is the value 
of 5vo at r = ra, is the radius of the Alfven critical point 
at which the solar-wind outflow speed U equals v'a, VAa is the 
Alfven speed at r = r^, and 



\\{r) 



n{r) 



(7) 



The above results can be applied to coronal holes and the 
fast solar wind with the use of the following model profiles 
for n, Bq, and U: 



3.23x10^ 2.51x10^ 1.85x10^ 



J5.6 



v.3.76 



cm 



-3 



Bo 



1.5 

2" 



Gauss, 



and 



U ^9.25x 10'2-cms"\ 
n 



(8) 



(9) 



(10) 



^ at 1 AU. Equations (O through ( fTOl l 
1.6O/?0, and t/(l AU) = 750 km s^^. 



where;*; = k/Rq, B is Bo in Gauss, and « is « in units of cm ^ . 
The density in equation ([8]) is the value from equation (4) of 
iFeldman et alj ( 119971) plus an additional r^^ component cho 
sen to give « = 4 cm 
givera = ll.liJ©, rm 

and lead to the U and v'a profiles shown in figure [T] The value 
of 5vo from equation (|5]l is also plotted in figure [T] where the 
value Sv'a = 155 km/s has been chosen so that 5vo remains 
bounded by the range of non-thermal line widths obtained by 
lEsser et aTl (Il999h multiplied by V2 to convert from an rms 
line-of-sight velocity to an rms velocity in the plane perpen- 
dicular to B. 



In Alfven-wave turbulence, fluctuation energy cascades 



from A- 1 ^ L I to smaller X i 



and then dissipates at A,x ^ Pp 



due to some combination of ion and electron damping. It is 
assumed that dissipation is negligible at A-x > 2pp, and that 



(11) 



at 2pp < A-x < ix, where hv^^ is the rms amplitude of the 
fluctuating velocity at perpendicular scale Ax. The con- 
stant a is related to the velocity power spectrum P^{k\^. 

If Pv °= kl"" for < fcx < (2pp)"', then a = (a - l)/2. 
Observations of sol ar-wind velocity flu ctuations at 1 AU 
find that o = 3/2 jPodesta et all 120071) . Numerical sim- 
ulations of magnetohydrodynamic (MHD) turbulence gen- 
erally find that g = 5/3 when 5v o ^ va fCho & Vishnia^ 
120001: iMiiller & BiskampI |2000; Haugen et al. 2004 ) and 
g = 3/2 when 5vo < 0.2va ( Maron & Goldreich 2001; 
Muller & GrapDin 2005; Boldvrev 2006; Mason et al. 200l 
iPerez & Boldvrev.2008. ,2009) . From figure[Tl 8vo/va < 0.3 
at r < 15^0, indicating that this near-Sun region is better 
described by simulations with 5vo < 0.2v'a than simulations 
with Sv'o ^ V'A. Based on the above studies, it is assumed 
that g = 3/2 and 

a ^ 0.25 (12) 

at r < 15Rq, the region on which this study focuses. 

UVCS observations show that has a larger thermal 
speed than protons at r ^ 2/?©, and hence a gyroradius that 
is several times larger than pp. Based in part on this obser- 
vation, it is assumed that pi > 2pp for minor ions and alpha 
particles at r > 2Rq. Equation ( flll l then gives 



5vi = ttiSvo 



Pi_ 



with 



ai = 1 (for He++ and minor ions). 



(13) 



(14) 



In the case of protons, electron Landau damping and stochas- 
tic proton heating drain energy from the cascade at Ax ~ pp, 
reducing ?>\\^ at Ax ^ Pp relative to the scaling in equa- 
tion (fTTT i. To account for this, it is assumed that 



ai ~ 0.71 (for protons), 



(15) 



where the particular value in equation ( fTSl ) is chosen so that 
the results in section |4] match observations of 7xp in coronal 
holes, as described further in the discussion of figure[3] Equa- 
tions ([B and ( HjI ) imply that 



Ei = ai 



6vo 

vxi 



l-a 



Lxiii 



(16) 



Equation ( fTSI l shows that as Txi decreases, Ei increases, which 
in turn increases the stochastic heating rate. One way of un- 
derstanding this is that Alfven-wave/KAW fluctuations with 
Ax ^ Pi cause the electrostatic potential <I> to vary in a compli- 
cated way in the plane perpendicular to Bo, with an rms varia- 
tion of 8<I>i ^ PiS/ii over a distance Pi, where hEi ^ Sv'iBo /c is 
the rms amplitude of the electric-field fluctuation at Ax = Pi- 
(The larger but smoother spatial variations in <I> associated 
with Alfven-wave fluctuations at Ax '3> Pi are ignored here, as 
they lead to drift motion rather than stochastic orbits.) Equa- 
tion ( fT3] l then gives 54>i <^ v^lI", and 5<I>i/v^j oc y'Jr^ For 
a < 1, decreasing vxi causes ^iSOi to become an increasingly 
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large fraction of the ions' perpendicular kinetic energy, and 
the ions' motion in the plane perpendicular to Bq becomes 
increasingly chaotic as a result. 

A possible objection to setting pi > 2pp at r > 2Rq for 
He++ and minor ions is that preferential heating of heavy ions 
is assumed from the outset, leaving open the question of how 
such ions first reach temperatures exceeding T±p. This initial 
evolution can be understood by considering He++, 0+^, and 
Fe+'^ ions with temperatures ^ T±p at small r. Given their 
charge-to-mass ratios, these ions have gyroradii that are com- 
parable to Pp when Txi ^ 7xp, and hence values of 5vi that 
are comparable to the proton value. As a result, the slower, 
heavier ions have larger Ei and much larger stochastic heating 
rates than protons when T±i ^ T±p. These larger heating rates 
then lead to T±i ^ T±p at larger r. 

In th e Cluster measurements at 1 AU analyzed by lBale et al.l 
( 12005b . the electric-field power spectrum is slightly larger 
at k±Pp = I than one would expect from an extrapolation 
of the power-law scaling that is present at smaller k±. The 
electric-field spectrum is a good proxy for the spectrum 
of (electron ) velocity fluctuations asso ciated with AW/KAW 
turbulence (ISchekochihin et all 120091) . and thus ai > 1 for 
protons in this data set. Nevertheless, the value tti = 
0.71 is reasonable for protons in coronal holes for two 
reasons. First, AW/KAW turbulence is more "imbal- 
anced" in coronal holes than at 1 AU, with a greater ex- 
cess of outward-propagating waves over inward propagating 
wave s (ICranmer & van BaUegooijenI l200l IVerdini & Vellil 
l2007h . Such imbalance weakens the cascade rate, causing the 
power spectrum to decrease more rapidly with increasing k± 
near k±Pp~l in response to a fixed KAW damping rate. Sec- 
ond, the linear damping rate of KAWs at k±Pp — 1 is signif- 
icantly larger in the P ^ 1 conditions of coronal holes than 
in the typi cal 3 ^ 0.5 — 1 condi tions found at 1 AU (see, e.g., 
figure 2 of lHowes et all (l2008ah '). 



4. ION TEMPERATURE PROFILES 

As described in the introduction, linear da mping of KAW s 
by ions is extremely weak when P <C 1 (lOuataertj 119981) . 
Based in part on this finding, it is assumed in this section that 
stochastic heating is the primary ion heating mechanism in 
coronal holes. The T±i profiles that result from this assump- 
tion are then calculated, and it is shown that for plausible pa- 
rameter values the resulting profiles provide a good fit to the 
observations. 

At r > 2Rq, collisional energy exchange between particle 
species can be neglected to a reasonable approximation, be- 
cause the e nergy exchange time scale exceeds the expansion 
time scale (lEsser et al.lll999h . 

fexp = ^. (17) 

The time scale fcond on which ion temperatures are modified 
by ion thermal conduction is > '"/v||i, where this lower limit 
corresponds to ions streaming freely at their parallel thermal 

veloci ty V||i = ^fa7 ||i/mi. At r > 2Rq, ion flows are super- 
sonic dKohl. J., et al.|[T998h . and thus fcond > fexp even in the 
free-streaming limit. Ion thermal conduction can thus be ne- 
glected at r > 2Rq to a first approximation. 
Because collisions and conduction are weak, minor ions at 



r > 2Rq evolve towards a state in which 

fh~fexp (18) 

for the following reasons. Minor ions draw very little power 
from the turbulence, and the power spectrum of the turbulence 
is essentially independent of the minor-ion heating rate Qi. 
Equations ([l]) and (|2]i thus imply that Qi is a highly sensi- 
tive function of the minor-ion temperature, with strong heat- 
ing at small T±i and exponentially weak heating at sufficiently 
large T±i. If fh ^> ?exp at some radius ri, then adiabatic cool- 
ing would cause T±i to decrease in proportion to Bq, and 
c/lnTxi/c/lnr would be < —2 near ri. The radial decrease in 
T±i near ri would lead to a sharp decrease in t^, and at larger r 
the minor ions would leave the adiabatic regime. Conversely, 
if fh <C ?exp at some radius r2, then dlnT±i/dlnr S> 1 near r2. 
The radial increase in Txi would cause fh to increase rapidly 
with increasing r, and at larger r the ions would again ap- 
proach a state in which fh ^ fexp- 

The observed 7xi profile of ions in coronal holes (see 
figure [3] below) appears to provide an example of the sec- 
ond case just described, in which ions swiftly evolve from 
a state in which fh ^ fexp to a state in which fh ^ fexp- The 
temperature is < 10^ K at r < 1.6Rq, presumably due 
to collisions with protons. As ions flow outward past 
r = 1.6R(v,, their temperature increases rapidly to ^ 10** K at 
r = 1.9Rq, where collisional energy exchange with protons 
is weak. At r = 1.9Rq, the 0+^ temperature profile abruptly 
flattens, with T±i remaining fairly constant out to 2.7/?0, indi- 
cating that fh ~ fexp at r > 2Rq. These rapid radial variations 
in T±i and dT^ Jdr present a challenge for theoretical models, 
but can be naturally explained in terms of stochastic heating 
— stochastic heating of an initially cool minor-ion population 
quickly increases 7xi, but then saturates at large T±i because 
of the reduction in orbit stochasticityQ 

The strong dependence of fh on E; implies that fh ^ fexp only 
within a narrow interval of £i values. The midpoint of this 
interval can be found by equating fh and fexp, which yields 

The right-hand side of this equation is ^ 1 in coronal holes, 
leading to values of £i that are ^ 1. For example, £l[r/U\ = 
6.13 X 10^ for ions at r = 2Rq given the assumptions 
listed in the caption of figure|2] Equation ( fT9l ) then gives Ei = 
2.96 X 10^^, assuming C2 = 0.15. For such small values of 
the value of Ei becomes relatively insensitive to changes in 
the right-hand side of equation ( fl9] l. For example, if Q.ir/Ui 
is increased from 6.13 x 10^ by 50%, the resulting decrease 
in Ei is only 5%. 

For protons, the comparatively flat T^i profiles seen in 
UVCS observations (Kohl et al. 1998) rule out the possibil- 
ity that fh ^ fexp or fh ^ fexp, assuming stochastic heating is 
the dominant heating mechanism. Thus, fh ^ fexp and pro- 
tons approximately satisfy equation ( fT9] l. However, protons 
can attain the required value of Ei not only by getting hotter 
or cooler, but also by absorbing energy from the turbulence 
and reducing the value of bvx^ at X± ^ Pp. No attempt is 

' It should be noted that the observed 0+' temperature profile has 
also been approximately reproduced by models invoking resonant cy- 
clotron heating by high-frequency Alfven/ion-cyclotron waves — see, e.g., 
Cranmer et al. 1 1999) and Isenberg & Vasquez (2009). 
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Fig. 2. — The values of Ej from equation 119) for protons, alpha particles, 
and O^^ ions, where ci = 0. 15 and Q.[ (r) is calculated using equation jg). Fo r 
protons, {/; in equation )19t is set equal to the value of U in equation jlOt . 
For He++ and U, is set equal to 1.75(7. 



made in this paper to treat the coupled evolution of protons 
and gyro-scale KAW fluctuations self-consistently. Instead, 
proton damping (and electron Landau damping) of fluctua- 
tions at X± ^ Pp is modeled simplistically by setting Ui = 0.71 
for protons in equation ( fTSl l. This particular value is chosen to 
match the UVCS observations shown in figure[3] as described 
further below. Although h elium comprises o nly ^ 5% of the 
ions in the fast solar wind (iBame et alJll977h . alpha particles 
are hotter than protons in the fast wind, and ma y also drain a 
significant amount of power from the cascade (Marsc h et al.l 
[1982; Kasper et al. 2007)0 The back reaction of helium heat- 
ing upon the turbulence, however, is neglected in this paper 

To determine £; from equation ( fT9] l. Bo is taken from equa- 
tion (|9]l and Ui for protons is set equal to U in equation (llOl l. 
For other ion species, Ui is taken to be 1.75 times the proton 
speed, which is consistent with UVCS observations of pro- 
tons and ions at r = 3Rp, (Kohl et al. 1998), but is only a 
rough approximation for other ion species and at other radii. 
The resulting values of ei for protons, alpha particles, and 0+^ 
ions are shown in figure|2]for 1.8/?0 < r < 157?0. Since ion 
thermal conduction and collisional energy exchange between 
particle species are neglected, the curves are not extended to 
smaller r where these processes become important. 

Once 8i(r) is determined using equation (T% . T±i can be 
determined from equation (fTSI l. which can be re-written as 



2k^ 



2/(1- 



As illustrated in figure |2] equation ( fT9] l leads to similar values 
of Ei for alpha particles and minor ions. If Ei and ai are the 
same for some set of ion species (or if tti/Si is the same), then 
equation (|20] | implies that 



'±i ' 



for these ions, where A — nii/nip, Z 



(21) 

qi/e, e is the proton 



^ For example, the Helium heating rate is comparable to the proton heat- 
ing rate for the temperature profiles in the right panel of figure [5] given the 
assumed values of U\ . 



charge, and I 
pressed as 



2fl/(l — a). The value of I can also be ex- 



l = 



2a-2 



(22) 



where o is the spectral index of the velocity power spectrum 
defined in section |2] In this paper, it is assumed that a — \/A 
and a = 3/2, which gives / — 2/3. If instead a = 1/3 and 
a = 5/3, then Z = 1. 

The largest radius at which equations ( fT9l ) and ( l20t can be 
applied is determined by the condition that the parallel ther- 
mal speed V||i remain ^ va so that ion Landau damping and 
transit-time damping of KAWs can be neglected. Assuming 
that «o and Bq are given by equations ^ and (|9]l, that 7||i('") 
for protons is no greater than the r-dependent (isotropic) 
temperature in equation (47) of ICranmer & van BallegooijenI 
(2005), and that protons are the ion species with the largest 
value of V||i, one finds that V||i < va/3 for all ion species at to 
at least 15Rq. 

Temperature profiles from equations ( fT9] l and ( |20l i are plot- 
ted in figure |3] The left panel of this figure shows Txi for 
protons and ions at 1.%Rq < r < 3.3Rq, and the right 
panel shows the T±i profiles for four ion species out to 15^0. 
The data points in bot h panels represen t the ion kinetic tem- 
peratures obtained by lEsser et alj ([1999) from UVCS obser- 
vations. The model temperature profiles in this figure were 
constructed using the «, Bq, and U profiles in equations ([HJ 
through ( fTOl ) and the values of L± and 5vo given in equa- 
tions ([H and ^ with 5v'a = 155 km/s as in figure [T] The 
value C2 = 0.15 was chosen to match the observed tem- 
peratures at r > 1.8/?0. The value ttp = 0.71 was then chosen 
to match the observed proton temperatures at r > 1 .^Rq. 

Possible sources of error in figure[3]include uncertainties in 
the values of C2, a, Ui, and t/i(r). For example, if C2 is in- 
creased from 0. 15 to 0.34, but all other model parameters are 
held fixed, then Txp decreases by a factor of 4.6 at r = 2Rq. 
If a is increased from 1/4 to 1/3 but all other parameters are 
fixed at their original values, then Txp decreases by a factor 
of 27 at r = 2Rq. If tti is increased from 0.71 to 1 for pro- 
tons, then 7xp increases by a factor of 2.5 at r = 2Rq. The 
value of Ui (r) has less of an impact on the ion temperatures 
that follow from equations (T% and ( |20| ). If U is decreased 
from \ .15U to \ .5U for 0+^ ions, then Txi increases by only 
5.2% at r = 27^0. 

There is, however, a potentially larger source of error as- 
sociated with Ui{r). As the relative velocity MJi between 
heavy ions and protons approaches va, the electric field in 
the heavy-ion frame decreases, since most of the Alfven- 
wave fluctuations propagate away from the Sun in the pro- 
ton frame. As a res ult, the s tochastic heating rate Qi for 
heavy ions decreases dChandra n et al. 2010). The accelera- 
tion of heavy ions to relative flow speeds ^ va provides a 
second possible mechanism (in addition to the reduction of £; 
through the increase of Txi) for saturating stochastic heat- 
ing, provided the amplitudes of sunward-propagating waves 
are much less than the amplitudes of anti-Sunward waves at 
A-x ~' Pi. The effects of MJi on 2±i are neglected in figure [3] 
which is justifie d at r < 3R p,, where UVCS observations show 
that AUi < Va dKohl. J., et al. 1998). However, for the veloc- 
ity profiles assumed in the construction of figure [3] AUi ^ ^a 
near r = 15Rq, suggesting that the heavy ion temperatures 
near r = 15Rq may be over-estimated by equation ( l20b . 

A self-consistency check on the assumed value of tti for 
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Fi g 3. — Left panel: the solid (dotted) line is 7^; for protons (O^^) from eq uati ons (19) and J20t . Right panel: the different curves show Txi from equations 1191 
and j20) for four ion species over a larg er ra nge in r. Both panels: equations )19t and \20\ are evaluated with a.\ = 0.71 for protons and Kj = 1 for other ions. The 
value of U\ is set equal to U in equation jlOt for protons and 1 .15U for other ions. For all curves, q =0.15 and a = 0.25, L± and 5i'o are taken from equations j4) 
and (5) with 5i'a = 15 5 km/s, and n and B are taken from equations (8) and (5). The circles and triangles correspond to observed kinetic temperatures for O^' 
ions and protons from lEsser et alj jl999l) . with the en'or bars included in the left panel only. 



protons can be obtained by comparing the cascade power T 
with the heating rate 2p,crit = ^B^xp/ ('Wp^exp) required to sus- 
tain protons at the temperature T±p given by equations iT% 
and ( l20l i. If the condition Qp.crit < T were not satisfied, then 
protons would absorb so much energy from the cascade that 
ttp would drop below the assumed value of 0.71, thereby re- 
ducing Qp.crit below r as required by energy conservation. 
The value of Qp.crit/r is shown in figure |4] for 1.87?0 < r < 
15Rq for the same set of assumptions used to obtain the 
proton temperature profile in figure |3] It can be seen that 
Gp.crit < r throughout this range of r. However, some caution 
is warranted here. As noted by Chandran & Hollweg (200^, 
equation (|6]l likely overestimates the turbulent dissipation rate 
for two reasons. First, large-scale Alfven waves launched 
from the Sun must propagate some distance into the corona 
before their energy cascades all the way to the dissipation 
scale. Second, equation ^ is derived in the limit of sm all L \ . 
As shown in figure 3 of Chandran & HoUwed (120091) (and in 
figure 5 of Dmitruk et al. (2002) for the case U = 0), finite 
values of L± reduce F relative to the small-Lx limit. These 
effects are particularly important close to the Sun. The small 
value of 2p.crit/r at r ^ 2Rq shown in figure|4]may thus sig- 
nificantly underestimate the fraction of the cascade power ab- 
sorbed by protons at this location. 

5. CONCLUSION 

Low-frequency Alfven-wave/KAW turbulence offers a 
promising explanation for the detailed features of ion kinetic 
temperatures seen in UVCS observations of coronal holes, 
including the abrupt radial variations in the tempera- 
ture profile, the widely different temperatures of protons and 
ions, and the large temperature anisotropy {Tj_[ 3> 
7]|i). When p ^ 1, ion heating from the linear damping of 
low-frequency Alfven waves and KAWs can be neglected. 



1 p 




0.01 I ' ' — ' — ' — I 

1 10 

r/R^ 

Fig. 4. — Op.crit is the approximate heating rate required to sustain protons 
at the te mper ature in the right panel of figure [5] given the outflow speed in 
equation ilOi . F is the turbulent dissipation rate in equation (6). 



and stochastic heating is arguably the primary way in which 
low-freque ncy Alfven-wa ve/KAW turbulence heats ions. As 
shown by Chandran et al.l fl2010 ). the stochastic heating rate 
is a strongly decreasing function of T±i. At small vxi, Ei is 
comparatively large and ion gyro orbits are strongly perturbed 
as the ions traverse the electrostatic potential of the gyro- 
scale fluctuations. This leads to chaotic ion orbits and strong 
stochastic heating. In contrast, at large v±i, ion gyromotion 
is only weakly perturbed by the electrostatic potential of the 
turbulent fluctuations, ions drift smoothly with nearly circular 
orbits in the plane perpendicular to B, and stochastic heat- 
ing is weak. As discussed in section |4] the observed 
temperature increases rapidly from ^ 10^ K to ^ 10^ K be- 
tween r = 1 .6Rq and r = 1 .9.^0, as the ions leave the region at 
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r < 1 .67^0 in which collisional energy exchange with protons 
is important dEsser et al.ll 19991) . At r = 1.9Rq, the tem- 
perature profile abruptly flattens, and Txi remains ^ 10^ K 
out to 2.1 Rq, the largest radius observed. This temperature 
profile is consistent with rapid stochastic heating of at 
7xi ^ 10^ K that saturates at 7xi ^ 10** K due to the reduction 
in orbit stochasticity. 

Ion thermal conduction, like collisional energy exchange 
between particle species, can be neglected to a good approx- 
imation at r > 2Rq. As a result, the strong temperature de- 
pendence of the heating rate causes all minor ion tempera- 
tures to evolve at r > 2Rqj to values at which fj, ^ fexpi at 
lower Txi, minor ions are rapidly heated, while at higher tem- 
peratures heating is ineffective and ions cool adiabatically. 
Because of the exponential dependence of fh on Ei, the con- 
dition fh ^ fexp is satisfied by different ion species at nearly 
equal values of Ei. This then leads to slightly more than mass- 
proportional temperatures — in particular, 7xi o<: A • (A/Z)' for 
ions with the same values of £; and tti, as in equations ( |20] | 
and (l2n i. The proton temperature is even smaller relative 
to minor-ion temperatures than the scaling in equation ( |20| | 
with constant tti would suggest, because proton heating and 
electron Landau damping reduce the amplitude of the fluctu- 
ations at A-x ^ Pp, which make the largest contribution to the 
stochastic proton heating rate. The condition fh ^ hx-p thus 
leads to temperatures that are much higher than the pro- 
ton temperature, consistent with UVCS observations. As de- 
scribed briefly in section |2] stochastic heating primarily in- 
creases the speed of ion thermal motions perpendicular to B 
when P <C 1 (Chand ran et al.ll2010i) . Because of this, stochas- 
tic heating also offers an explanation for the observation that 
rxi > for io ns (iKohl. J., et al.lfT998t iLiet al.lfT998t 
lAntonu cci et al. '2000V 

In section |4] 7xi profiles are calculated for several ion 
species from the condition ~ fexp using an observation- 
ally constrained model of Alfven-wave turbulence in coronal 
holes. For plausible values of the model parameters, the re- 
sulting temperature profiles provide a good match to obser- 



vations of protons and ions. However, there are sev- 
eral sources of uncertainty in these calculations, the largest 
of which are associated with the parameters C2, a, and ttp 
(the value of tti for protons). The constant C2 relates to 
the efficiency of stochastic ion heating and was evaluated by 
Chandran et al. (2010) using numerical simulations of test- 
particle protons interacting with randomly phased waves in a 
low-p plasma. However, the value of C2 for strong turbulence 
and for other ion species is not yet known. The quantities a 
and ttp describe, respectively, the slope of the inertial-range 
velocity power spectrum and the decrement in the velocity 
spectrum at A-x ^ Pp arising from electron and ion damping. 
In order to determine whether stochastic heating can indeed 
explain the ion temperatures observed in coronal holes, future 
work is needed to determine C2, a, and ttp more accurately. 
For exampl e, direct numerical simulations or a cascade model 
(see, e.g., Howes et all (12008a)) of strong AW/KAW turbu- 
lence including both stochastic heating and electron Landau 
damping would help ascertain the value of ttp. Direct nu- 
merical simulations of strong, highly anisotropic, AW/KAW 
turbulence interacting with test particles could be used to de- 
termine the value of C2 that is appropriate for coronal holes. 
In addition, in situ measurements from NASA's planned Solar 
Probe mission will establish the power spectrum of AW/KAW 
turbulence and determine Ui{r) and 7xi(r) for several ion 
species at heliocentric distances as small as ~ 9.5Rq. These 
measurements will lead to a rigorous test of the stochastic 
heating model, including the 7xi(r) predictions shown in fig- 
ure[3]and the dependence of 7xi on ion mass and charge given 
in equation d^TT i. 
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